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Abstract 

We consider a system of ultracold atoms in an optical lattice as a quantum 
simulator for electron-positron pair production in quantum electrodynamics 
(QED). For a setup in one spatial dimension, we investigate the nonequilib¬ 
rium phenomenon of pair production including the backreaction leading to 
plasma oscillations. Unlike previous investigations on quantum link models, 
we focus on the infinite-dimensional Hilbert space of QED and show that 
it may be well approximated by experiments employing Bose-Einstein con¬ 
densates interacting with fermionic atoms. Numerical calculations based on 
functional integral techniques give a unique access to the physical param¬ 
eters required to realize QED phenomena in a cold atom experiment. In 
particular, we use our approach to consider quantum link models in a yet 
unexplored parameter regime and give bounds for their ability to capture 
essential features of the physics. The results suggest a paradigmatic change 
towards realizations using coherent many-body states for quantum simula¬ 
tions of high-energy particle physics phenomena. 

Keywords: 


Introduction. The creation of electron-positron pairs from the vacuum 
of quantum electrodynamics in an external electric held is a longstanding 
prediction that has not yet been directly observed meg. Upcoming ex- 
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perimental laser facilities, such as the Extreme Light Infrastructure (ELI) 
[3j, start approaching the required critical field strength of E c ~ 10 16 V/cm. 
Theoretically, the non-linear interplay of the produced many-body states and 
the applied field represents a remarkable challenge with important links to 
a wide range of elusive phenomena such as Unruh and Hawking radiation or 
string-breaking in quantum chromodynamics (QCD) [3HE]. 

Whilst the critical field strength E c = M 2 /e is determined by the elec¬ 
tron/positron mass M and the absolute value of the electric charge e in 
QED, the pair-production phenomenon essentially depends on the dimen¬ 
sionless ratio E/E c > 1 for an applied electric field E. In principle, physical 
systems with very different characteristic scales can thus be used to realize 
the underlying phenomenon. It has recently been suggested to employ a 
system of ultracold atoms in an optical lattice to study the physics of pair 
production and string breaking m Even though the implementation of a 
gauge symmetry in an atomic setup is demanding in general pm, it may 
provide a unique way of answering crucial open questions, such as regard¬ 
ing the nonequilibrium dynamics of the strong and electroweak sector of the 
standard model of particle physics probed in heavy-ion collision experiments 
or early-universe cosmology H2 US- 

Many proposals concentrate on quantum link models [33] rather than 
QED, or even QCD. Since the Hilbert space of a quantum link model is 
finite-dimensional, the mapping to atomic systems is expected to be greatly 
facilitated. However, it is a crucial question how much of the physics of 
the infinite-dimensional representation corresponding to QED may be cap¬ 
tured in practice. Theoretical estimates based on diagonalization or matrix 
product states techniques are typically limited to low-dimensional represen¬ 
tations HSHH]. Recently, powerful functional integral (FI) techniques have 
been employed to simulate the real-time dynamics of pair production and 
string-breaking directly in QED on a one-dimensional lattice pT8Ll2(Tj . and in 
three dimensions [21] • This progress has become possible since strong bosonic 
fields can be efficiently sampled from coherent classical fields while keeping 
the full quantum nature of fermions [ T8I [22 j. 

In this work, we exploit this observation and start from the infinite- 
dimensional representation of the QED gauge group, pointing out that it 
may be well approximated by experiments using Bose-Einstein condensates 
interacting with fermionic atoms. For the first time, by using FI techniques 
we can estimate the physical parameters required to describe the QED phe¬ 
nomenon of pair production in a cold atom setup and we discuss the exper- 
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imental realization. For this paradigmatic example, numerical studies using 
the FI approach are still feasible and serve as an important benchmark for 
future quantum simulation experiments in numerically inaccessible regimes. 
In particular, we use our approach to consider quantum link models in a yet 
unexplored parameter range and to give bounds for the dimensionality of the 
employed representation in order to capture essential features of strong-field 
QED in those formulations | jl6] . 

While the FI techniques can also be applied in higher dimensions and 
non-Abelian gauge theories |21j . we focus here on the conceptually important 
example of QED in one spatial dimension |23, 23J- Since there are no spatial 
plaquette terms in this case, angular momentum conserving atomic scattering 
processes can be used to directly implement the [/(1) gauge symmetry p3j. 
From a phenomenological point of view, this theory shares key properties 
with QCD as, e.g., dynamical string breaking, and hence provides valuable 
insights into nonequilibrium aspects of the theory of strong interactions. 

Cold atom gauge theory. We start with the Hamiltonian formulation of 
lattice QED using the staggered fermion discretization m- To this end, 
the spinors are decomposed such that particle and antiparticle components 
separately reside on two neighboring sites of the lattice. Introducing the 
operators for the fermion held i/j n , the link U n being connected to the gauge 
potential A n and the electric held E n , the Hamiltonian reads 

^QED = + 

n 

~ Y a ^nUn^n+l ~ V’l+i U l^n 

where a is the lattice spacing and g the gauge coupling. The dynamical 
variables of QED fulfill [E n , U m \ = gS nm U m . The Gauss-law operator G n = 
E n — E n _ i — gij^n commutes with the Hamiltonian [G ni -Hqed] = 0. The 
last relation manifests local gauge invariance. 

To make contact with quantum link models, we approximate QED by 
another theory, which is still gauge invariant and allows for an experimental 
realization with cold atoms We substitute E n —> gL z , U n — > [l(£ + 

l)] _1/2 L +iri , where (L x>n , L y ^ L ZjU ) are quantum spin operators which obey 
[L i)7l , Lj.m] = i5 nm e ijk L ktn , and the raising operator is L +}U = L x , n + iL y<n . 
The performed substitution renders the dimension of the local Hilbert space 
hnite-dimensional. Consequently, the QED relation [U n , J = 0 is no longer 
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valid but is replaced by [L +>n , L_ >m ] = 2 8 nin L Ztin . Remarkably, local gauge 
invariance is not affected. 

The Schwinger representation of the angular momentum operators allows 
us to express the quantum spins by bosonic degrees of freedom b n and d n J26]: 
L +)U = b^dn, = d^ n b n and L ZjU = — d) n d n )/ 2, with the constraint 

fof b n + d) n d n = 2 L Here, £ denotes the spin magnitude fixing the number of 
bosonic atoms. In this representation, the Hamiltonian describing the cold 
atom (CA) system becomes 

hca = 

} ( 2 ) 

and depends on two species of bosonic operators b n , d n living on links and 
fermionic operators ^ n located on lattice sites. The parameters are deter¬ 
mined by the basic physical quantities of the cold atom system, such as the 
gauge coupling g given by the on-site scattering processes of the bosons and 
the spin-changing collisions between the fermionic and bosonic atoms. 

Employing the density-phase representation with b n = \J £ + 5pb :n e l6b ’ n 
and d n = £ + 8pd,nd l6d ' n , one finds H C \ = -Hq ED + 0(8p/£). Therefore, the 

number of bosonic atoms per site, 2£, controls the approximation and allows 
to tune from the quantum link formulation to lattice QED by increasing £. 

In previous work, the Hamiltonian (|2]) was studied for £ ~ 0(1) via 
diagonalization or matrix product states methods PUMHET]. Here we 
consider for the first time the dynamics in the regime £ 1 to approach 

lattice QED. From an experimental point of view, this regime corresponds 
to putting Bose-Einstein condensates on the links rather than single bosonic 
atoms. 

Functional integral approach. To study the strong-field regime of QED, 
the field strength needs to be of the order of the critical field E c = M 2 / g. 
The corresponding cold atom setup is characterized by E c = g \ N h — N d \ /2 ~ 
M 2 /g r, where Nf } , N c i denote the number of atoms in the Bose-Einstein con¬ 
densates. For Nb,Nd ~ 0(£) 1, the FI approach of Refs. 11.8.1201-1221 [28 j 

allows us to study the dynamics in this regime. 

To this end, we denote the bosonic fields collectively by (ft n = (b^, b n , d^ n , 
and define the generating functional for correlation functions in the presence 
of sources J n = {J b ,n, </£„, J d , n , J^n) by Z[J] = Tr{p 0 T c e jJ '^}. Here, p 0 is the 
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initial density matrix, J ■ 0 = Yin It Jn(t) ' <Pn{t) with t the time coordinate 
along the closed time path C, and Tq denotes time-ordering along C. Em¬ 
ploying the coherent state basis, the matrix element of the density operator 
at initial time is (+|po | — ), where |+) and |—) are the Erst coherent states 
on the forward and backward contour, respectively. The FI representation 
of the generating functional is 

Z[J\ = I [#][#*#] (+| p 0 |-> e iS+iJ ■* (3) 


with the action 

S= j 'Y^ n id t ‘ij) n + b ] n id t b n + d\id t d n) ~ H C A ■ 


( 4 ) 


We analytically perform the Gaussian integral for the fermions and then 
expand to first order in the bosonic response held <j) n , which arises from 
the Keldysh rotation 0 n = 0 n + sgn c (f> n . Disregarding higher-order terms, 
i.e. neglecting subleading corrections in bosonic occupancies [21] that are 


suppressed by N b 1 , Ab 1 -C 1, leads to the self-consistent set of equations 


■ n i 9 ® jf r 7 . • d n F n + \ n 

io t b n = - b'b n b n + % - . , 

2 Aay/£(£ + 1 ) 

■Q , _ 9^1 +77 • b n F nn .)_1 

‘ ~ ~ " “ \a^e(i+ 1) ’ 

id t F nm = ^ [hnn'Fn'm ~ F nn >h^] . 


( 5 ) 


Here, F nm = i s the fermion equal-time correlation function, whose 

evolution is governed by 

1 TSnm. 

The equations (|5]) preserve the Gauss law, d t G n = 0, if initialized accordingly. 
A similar derivation for the QED Hamiltonian ([!]) gives 

dtEn = Y a MFn + i n u n l 
dtU n igaE n U n , 

id t F nm = , (6) 

n' 


h CA = 

' nm 


n—l^n— 1 m IjdJ n +1 m\ 

2a^/£(£ + 1) 


+ M(- 
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time: t ■ M 

Figure 1: Time evolution of the fermion number density in the cold atom system for 
different numbers of atoms ~ £ as compared to the QED result (t —> oo). The straight line 
corresponds to the Schwinger formula, which neglects the backreaction of the produced 
fermions on the applied field. 


with 


hZ D = ^[K-A-im - uj n+lm } + M(—l) n 8 nm . 

In fact, by taking the time derivative of E n —* g(b' n b n — d' n d n )/2 and U n —* 
b^dn and inserting the density-phase representation, one finds that ([5]) ap¬ 
proximates (J6]) with a truncation error of 0(8p/i). In the following, we keep 
i as a parameter to answer the question whether essential properties of QED 
can be captured for finite l. 

Pair production. The creation of electron-positron pairs in a uniform elec¬ 
tric field E may be viewed as a quantum process in which virtual electron- 
positron dipoles can be separated to become real pairs once they gain the 
binding energy of twice the rest mass energy. This QED process has been esti¬ 
mated una neglecting the backreaction of the produced pairs on the applied 
field, and the analytic result for the rate n = M 2 E/(2 tiE c ) exp(—nE c /E) is 
depicted in Fig. [lj 

This estimate should be valid at sufficiently early times and provides an 
important benchmark for any simulation method. Therefore, we consider 
a spatial lattice of length Na with periodic boundary conditions and leave 
more refined estimates taking into account specific trap geometries of cold 
atom systems for further studies. We first compute the real-time evolution 
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Figure 2: Dynamics of the homogeneous electric field as represented by the bosonic species 
population imbalance of the atomic system for different values of i as in Fig.[l] The plasma 
oscillations arise from the backreaction of the produced pairs on the applied field. 

according to (J5]) for g/M = 0.1 and aM = 0.005 in the limit i — > oo, 
where it agrees to QED described by (|6]) . Here N determines the number of 
fermionic atoms, and we checked that for the largest employed lattices with 
N = 512 no significant volume dependence can be observed and our results 
are insensitive to changes in the lattice spacing. Accordingly, employing a 
standard definition of the particle number density PSl El] the simulation 
result for QED (£ = oo) as shown in Fig. [I] agrees well with the Schwinger 
formula at early times. At later times the backreaction of the produced pairs 
on the applied field is seen to give the expected sizable corrections [IS, IPj- 

Of course, in the corresponding cold atom system no particles are pro¬ 
duced as the number of atoms is fixed. However, since two neighboring 
fermions are considered as particles and antiparticles in the staggered formu¬ 
lation, pair production is encoded in the hopping of atoms between odd and 
even sites of the optical lattice. This generates correlations, whose time evo¬ 
lution describe the corresponding phenomenon of pair production as shown 
in Fig. [l] The results demonstrate the convergence of the atomic system’s 
dynamics to the QED behavior as the number of bosonic atoms is increased. 
For i = 2500 we still observe considerable deviations from the QED result, 
whereas the difference becomes small for i = 5000. 

In Fig. [2] the time evolution of the cold atom analogue of the electric 
held, E = g(Nb — Nj)/ 2, is given for different values of i. As for Fig. [lj we 
start with a bosonic species imbalance Nb~ Na = 2 M 2 /g 2 > 0 corresponding 
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Figure 3: Momentum distribution of the produced fermions at tM = 25 (dotted), tM = 75 
( dashed ) and tM = 175 (solid) for l = 10 4 such that the QED result is well described. 


to the critical electric field strength in QED, and the analogue of the Dirac 
vacuum or ‘Fermi sea’ with the lowest Nj 2 energy eigenstates occupied. By 
comparison to Fig. [lj we observe a decrease of the electric field as the fermion 
number increases due to pair production. The correlated hopping of fermions 
reduces the bosonic species imbalance until it becomes zero and even changes 
its sign with N b — N d < 0 giving rise to plasma oscillations [18], [21, 29]. At 
times when the corresponding electric field drops below a critical level, par¬ 
ticle creation effectively terminates as reflected in the characteristic plateaus 
in the particle number density. 

The plasma oscillations are caused by the production and subsequent 
acceleration of particle-antiparticle pairs, which can be observed from their 
momentum distribution H5H2D as shown in Fig. [3] for £ = 10 4 such that the 
QED result is well reproduced. The homogeneous electric field dominantly 
produces fermions around zero momentum, and accelerates the particles to 
higher momenta during the subsequent evolution. As the electric field de¬ 
creases due to energy conservation, the production amplitude around zero 
momentum drops as well. At a time tM ~ 175, the fermions reach their 
maximum momentum along with a vanishing net electric field. Subsequently, 
the fermionic current results in a further decrease of the electric field to neg¬ 
ative values along with a deceleration of the produced particles leading to 
the observed plasma oscillations. 

The limiting experimental resources enforce a study of the dependence 
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time: t ■ M 


Figure 4: Time evolution of the electric field for different total numbers of fermionic atoms, 
N = 32 (dotted), N = 128 ( dashed ) and N = 512 (solid) with £ = 10 4 . 

on the total number N of fermionic atoms. Fig. [4] shows the electric field 
for different N with fixed l = 10 4 . We observe a reasonable description of 
the QED results for a full oscillation period employing a total number of 
fermionic atoms down to about N = 128, with sizable deviations occurring 
at later times. For N = 512 accurate descriptions are achieved for the entire 
range of times we considered. 

Experimental realization. The physics of QED pair production in one 
spatial dimension may already be realized with available experimental tech¬ 
niques [81 HUSHES]. Here we point out how the relevant regime of large 
i may be efficiently implemented and manipulated experimentally with the 
help of coherent many-body states. 

Regarding the bosonic degrees of freedom, we confine two substates of 
one hyperhne manifold of bosonic atoms in a one-dimensional geometry via 
an external potential. Further, we employ a red detuned laser to generate 
an optical lattice such that the atoms are localized and the nearest neighbor 
hopping is suppressed. Already this construction allows us to realize meso¬ 
scopic bosonic gases with two components (b n and d n ) per site which can be 
described by the one axis twisting Hamiltonian [32] corresponding to the 
first two terms of (|2]). The preparation of the bosonic atoms mimicking the 
electric held can be achieved by a magnetic held or homogeneous two-photon 
microwave coupling, which produces a coherent spin state between the two 
components with a non-vanishing population difference. 
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Concerning the fermionic degrees of freedom, we again trap two substates 
of one hyperfine state manifold of fermionic atoms in a one-dimensional ge¬ 
ometry. The aforementioned optical lattice is blue detuned for the fermions. 
Consequently, the fermions are located between the bosonic links, which is 
an essential ingredient of (j2]). In addition, we superimpose a second optical 
lattice with double lattice spacing in order to generate the staggered struc¬ 
ture of the fermions. The frequency of this laser is tuned closer to resonance 
with respect to the fermions than to the bosons, such that the second lattice 
does not affect the latter ones. The staggered structure leads to a mini-gap 
in the dispersion relation of the fermions such that the initial state of the 
fermions corresponds to a fully filled lowest mini-band [33]. The detection of 
the fermions can be achieved by exploiting the band mapping technique [33]. 

The overlap of neighboring bosonic and fermionic atoms makes the hop¬ 
ping of one fermion from one site to the next possible via spin changing 
collisions. This modifies the internal state of the boson, i.e., species b be¬ 
come d or vice versa. This interaction induced hopping process implemented 
by boson-fermion spin-exchange collisions locally preserves the total spin |35j. 
In particular, the detuning of this process that is experimentally controlled 
by the external held corresponds to the mass term in (|2]) . The dynamics 
of the system is initiated with a quench of the mass term from being far 
off-resonant. 

Conclusion. Simulating high-energy physics by experiments with ultra¬ 
cold atoms may be achieved with coherent many-body states rather than 
single atoms. Our findings of the required resources in terms of atom num¬ 
bers and protocols may already be realized with available experimental tech¬ 
niques. This opens new possibilities to resolve questions in the strong¬ 
coupling regimes of gauge theories, where no alternative real-time simulation 
techniques are known so far. Together with recent experimental proofs of 
concept in strongly interacting systems [36], one may hope to realize the old 
dream of solving complex problems in quantum held theory by experiment. 
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